Self-refraction of supernova neutrinos: mixed spectra and three-flavor instabilities 
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Neutrinos in a core-collapse supernova undergo coherent flavor transformations in their own back- 
ground. We explore this phenomenon during the cooling stage of the explosion. Our three-flavor 
calculations reveal qualitatively new effects compared to a two-flavor analysis. These effects are es- 
pecially clearly seen for the inverted mass hierarchy: we find a different pattern of spectral "swaps" 
in the neutrino spectrum and a novel "mixed" spectrum for the antineutrinos. A brief discussion of 
the relevant physics is presented, including the instability of the two-fiavor evolution trajectory, the 
3-flavor pattern of spectral "swaps," and partial nonadiabaticity of the evolution. 
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1. Introduction and motivation. - Corc-collapse su- 
pernovae play an essential role in the evolution of the 
Universe, from controlling the temperature of the gas and 
the rate of star formation in the galactic disk {e.g., [T]), 
to synthesizing and dispersing heavy elements {e.g., [5]). 
The ashes of ancient explosions have literally shaped the 
world around us. The struggle to understand and model 
the explosion mechanism has occupied researchers for the 
better part of the 20th century [3lfl2] and remains a very 
active topic (for review, see, e.g., P^fTl] ). 

Neutrinos, which are emitted during the first ~ 10 sec 
after the onset of the collapse and arrive to us directly 
from the core of the star, could serve as a direct probe of 
the explosion mechanism. It is expected that, unlike su- 
pernova 1987a [151 116| , the next galactic supernova may 
yield O(IO^) antineutrino {D^) events and, furthermore, 
a similar number of neutrino (I'e) events, if a large liquid 
argon detector is constructed at DUSEL. With such high 
statistics, it should be possible to reconstruct the second- 
by-second evolution of the and D^, spectra. The task 
is to understand how to "read" this signal, i.e., how to 
extract signatures of various physical processes from it, 
and to optimize the detector design. 

One important process that needs to be thoroughly 
understood is the coherent flavor transformations of neu- 
trinos outside of the neutrinosphere. Compared to the 
well-studied case of solar neutrinos, in a supernova these 
transformations are much more involved. Not only neu- 
trinos and antineutrinos of all flavors are emitted, not 
only are there two mass splittings - "solar" and "atmo- 
spheric" - to worry about [T7], but the physics of the 
transformations is significantly richer. For example, sev- 
eral seconds after the onset of the explosion, the flavor 
conversion probability is affected by the expanding shock 
front [H] and the turbulent region behind it ^9j. The 
conversion process in such a "bumpy," stochastic profile 
is qualitatively different from the adiabatic MSW effect 
in the smooth, fixed density profile of the Sun. 

Even more complexity is brought about by the coher- 
ent scattering of neutrinos off each other [201 - I36] . This 
neutrino "self-refraction" results in highly nontrivial fla- 



vor transformations close to the neutrinosphere, 

typically within a few hundred kilometers from the cen- 
ter, where the density of streaming neutrinos is very high. 
Since the evolving flavor composition of the neutrino flux 
feeds back into the oscillation Hamiltonian, the prob- 
lem is nonlinear. Furthermore, as the interactions cou- 
ple neutrinos and antineutrinos of different flavors and 
energies, the oscillations are characterized by collective 
modes. This leads to very rich physics that has been the 
subject of intense interest over the last several years. 

One may wonder whether all this complexity will im- 
pede the extraction of useful information from the future 
signal. In fact, the opposite is true: the new effects can 
imprint information about the inner workings of the ex- 
plosion on the signal. For example, by observing the sig- 
natures of the expanding shock and the post-shock region 
in the neutrino signal, we will learn about the develop- 
ment of the explosion during the crucial first 10 seconds. 
This information could be inaccessible in other ways. 

On the other hand, it is fair to say that the subject 
is still far from being exhausted and qualitatively new 
effects continue to be uncovered. This letter is another 
contribution to this effort. We explore the self- refraction 
phenomenon for the conditions typically present during 
the cooling stage of the protoneutron star (several sec- 
onds into the explosion) and show that a three-flavor 
analysis of this process reveals several new effects. Below 
we report the main results; a detailed discussion will be 
given elsewhere [49] . 

2. Formulation of the problem. - The energy spec- 
tra of emitted neutrinos evolve as the explosion devel- 
ops. Since the self-refraction effect is nonlinear, differ- 
ent initial spectra may result in qualitatively different 
outcomes. Therefore, it is important to choose spec- 
tra that are (i) specific to the late phase and (ii) sim- 
ulated with sufficient accuracy. It has been shown [SO] 
that the emerging spectra of the nonelectron neutrinos 
{vx) are sensitive to several seemingly subdominant pro- 
cesses, such as i^xe^ scattering, neutrino coannihilation, 
bremsstrahlung, and nucleon recoil. The addition of nu- 
cleon recoil, in particular, makes the i>x spectra softer and 
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FIG. 1: Spectra of Ue (left) and z/e (right) sit 1000 km, shown by the filled regions. The top plots are computed within the 
two-flavor framework, the bottom ones with a full three-flavor calculation. Inverted mass hierarchy and S13 = 0.01 are assumed. 
Also shown, by dashed/dotted curves, are the original spectra at the neutrinosphere, as labeled. 



brighter, changing the average energy from 26 MeV to 
~ 19 MeV [SU] [S3]. This significant shift may aher the 
collective oscillations, and needs to be included. 

Indeed, collective oscillations with such spectra were 
recently studied in [IS] and several interesting features 
were observed. In particular, for the inverted hierarchy 
(IH), the neutrino spectra were found to be exchanged 
between v,, and Vx, but only in an interval of energies, 
from about 6 to 22 MeV. The neutrinos outside of this 
interval were unchanged, resulting in the presence of two 
spectral splits. While spectral splits had been observed 
with other fluxes, starting with the seminal work [39j . and 
in fact seem quite ubiquitous, multiple splits is a new and 
interesting phenomenon, with potentially important im- 
plications for signal detection (high-energy splits would 
be easy to detect). 

While the analysis in [48] is an important step toward 
understanding the collective oscillations of late-time neu- 
trinos, it was limited to two flavors. What happens with 
all three flavors included in the calculation? A priori, 
several possibilities come to mind. The effects of the 



third state could be a small correction, as in the cases of 
solar, KamLAND, and atmospheric neutrino oscillations. 
Alternatively, the effects of the two splittings could be 
large, but factorizable, just like the conventional MSW 
effects in a supernova, which can be treated "pairwise" 
between the states and combined at the last step |17| . 
Instead, something even more interesting happens: the 
three-flavor evolution with a nonzero solar mass splitting 
gives an entirely new result, as we display next. 

3. Setup of the calculation. - For the late-time spec- 
tra, we use the results of the Monte-Carlo simulations 
from [50]. For the calculations shown in this letter, we 
select the point p — 10, q = 3.5 from Table 6 in that 
paper, corresponding to the spectrum emitted from a 
steep power-law profile near the neutrinosphere (p oc r^^, 
T oc r~'). For the matter profile at r ^ 100—1000 km we 
assume a neutrino driven wind with p = po(10 km/r)'^. 
We take po = 2 x 10^ g/cm'^, and = 0.5. 

Our three-flavor calculation is carried out with the fol- 
lowing parameters: |Am^tml = 2.7 x 10~^ eV^, AtUq = 
7.7 X 10-5 eV2, = 0.01, and sin^ 6*12 = 0.31. In the 
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FIG. 2: Impact of decreasing the solar mass splitting, Attiq, 
on the neutrino spectra at 1000 km. Notice that, while a 
strictly vanishing ArriQ gives the two-flavor result, even a 
tiny nonzero value of Ajtiq qualitatively changes the answer. 



FIG. 3: Impact of increasing the solar mass splitting Attiq on 
the neutrino spectra at 500 km. Large Athq makes the evo- 
lution more adiabatic, leading to the formation of a (broad) 
spectral split. 



two-flavor calculation, we set the solar mixing angle 9i2 
to zero and drop the state that in vacuum is separated 
from the predominately (i^e) state by Attiq. 

We perform multienergy, single-angle calculations of 
the evolution, starting at 40 km and ending at 1000 km. 

4. Results: comparison of two- and three-flavor runs. 
- The resulting spectra at 1000 km for the IH case 
(Am^jjjj < 0) are presented in Fig. [l] We concentrate 
here on this case, because it most clearly illustrates the 
effect; both hierarchies are treated and contrasted in [i^ . 
The top panels show the two-flavor calculations, the bot- 
tom ones, the corresponding three-flavor runs. The 
spectra are on the left, those for are on the right. 
The dashed curves and the dotted curves show the corre- 
sponding initial spectra (see legend). Animations show- 
ing the evolution of the spectra as a function of the radius 
for both IH and NH are available as Supplemental Ma- 
terial online (STJ [55] . 

The results of the two-flavor calculations appear to be 
in very good agreement with the inverted hierarchy cal- 
culations of |3H]. The three-flavor calculation results are 
visibly different: (i) the high-energy split in the Ve chan- 
nel is gone; (ii) in the Vf, channel, the flavor swap prob- 
ability is neither zero, nor one, but increases gradually 
with neutrino energy. The swap is incomplete in this 
case. 

5. Discussion. - Both of these results deserve further 
investigation. First we can establish that the physics 
driving flavor conversion before 1000 km are the neutrino- 
neutrino interactions, with the conventional MSW effect 
being unimportant. The atmospheric level crossing does 
occur here, but for the chosen parameters it is strongly 
nonadiabatic (flavor preserving). Moreover, it occurs 
when r > 600 km, by which point the neutrino self- 
refraction effects have ceased. The small MSW effects 
are seen in the 9^ channel as small wiggles. 



It looks paradoxical at first that the solar splitting 
completely changes the evolution of neutrinos at high en- 
ergies, even though it is only ~ 3% of the atmospheric 
splitting. Clearly, at this level Atoq cannot "overpower" 
Am^jjj^; what it can do, however, is change the nature 
of the collective motion, allowing it to extend from the 
2-flavor subspace into the full 3-flavor space. Evidently, 
this results in a radical change, not a small correction. 

Let us investigate how this change happens as the value 
of the solar splitting is turned on. The results are shown 
in Fig. [2] These at first may be even more surprising: the 
two-flavor spectrum is reproduced only when Atuq = 0; 
as soon as it is nonzero, even very small, the high-energy 
split disappears. Since for Atoq — 7.7 x 10"'' eV^ (1% 
of its true value) the corresponding oscillation length is 
10* km - much longer than the scales in the problem - 
one might think the two-flavor limit should be reached. 
Instead, the spectrum in this case is closer to the realistic 
three-flavor one than to the two-flavor one. 

To understand this result, it proves useful to examine 
the evolution as a function of radius in the matter mass 
basis [49l EH [52] • Neutrinos, initially in the mass eigen- 
states, develop an instability which leads to large collec- 
tive oscillations [571 133 [Mj- What is interesting in our 
case is that, shortly after the collective oscillations de- 
velop between the "atmospheric" eigenstates {1^2,3), the 
third state (I'l) joins in. Therefore, not only is the initial 
configuration unstable, but the two-flavor trajectory is 
also unstable in the three-flavor space. A small nonzero 
Attiq is enough to displace the system from the two- 
flavor subspace; the instability makes it then run away 
from it (driven primarily by Am^^^^). The role of small 
Attiq is thus similar to that of small 6*13 in the develop- 
ment of the Kostelecky-Samuel instability [271 [571 [55] . 

Since the evolution is no longer constrained to the two- 
flavor subspace, its final state may be expected to be 
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different. Indeed, it is. In the two-flavor case, the final 
state has a swap V2 -(-t- on the interval 6 < -B^ < 
22 MeV[54j. In contrast, the final state of the three- 
flavor calculation is significantly more complicated, as 
illustrated in Figs. |4]and[5] The pattern is the same 
only up to El, ^ 10 MeV, with the third state vi being 
a spectator, vi ^ vi. For 10 < £',y < 20 MeV we instead 
observe a cyclical permutation: — > J^s, I's — >■ i^i, J^i — >■ 
V2\ while for E^, > 20 MeV it is the ly^ state that is a 
spectator and the other two swap spectra: vi -H- 1^2, 1^3 ^ 
1^3. Since the initial spectra of ^1 and are identical, no 
second split is seen at 1000 km. See [H] for more details. 

Figure [2] also shows another important role of Atoq : 
as it nears its physical value, the spectrum at high energy 
becomes closer to the emitted Vx spectrum. This means 
the neutrinos at the end of the collective oscillations are 
"put" into the Hamiltonian eigenstates. This happens 
when the evolution is adiabatic [321 IH] HI]. Observe 
that adiabaticity is broken even when Am^ is 20-50% 
of its physical values. This means that adiabaticity for 
physical Attiq is only marginal in this channel. 

Indeed, let us compare the neutrino vacuum oscillation 
length to the scale height of the neutrino-neutrino poten- 
tial. The latter is a power law oc r^"*, so the scale height 
is \dliiH„„/dr\-'^ ~r/4~ 75-100 kmforr- 300-400 
km. The atmospheric splitting for E ^ 15 MeV gives a 
characteristic scale 2E / Am^^^ ~ 2 km, so a high degree 
of adiabaticity is expected (and seen for the low-energy i^e 
split). In contrast, for the solar splitting, 2E/ Arn^ ~ 77 
km, so the evolution is only borderline adiabatic |56| . 

The weakness of adiabaticity also helps to understand 
the Ve spectra: the antineutrinos are not placed into mass 
eigenstates. To increase adiabaticity, one can artificially 



increase the value of Ar 



We reran the calculation 



with Awq 5 and 10 times larger than its actual value. 
The resulting spectra, shown in Fig. [3j indeed exhibited 
a more "conventional" split, centered around E^ ^ 19 
MeV, although still fairly broad 5 MeV and 3 
MeV half- widths correspondingly). The mixed spec- 
trum found for physical Attiq can be thought of as an 
extremely broad split, with a width comparable to the 
entire range of the antineutrino energies. (The width of 
the split is related to the degree of adiabaticity, [H].) 

6. Conclusions. - The late-time spectra provide an in- 
teresting physical system for studying collective transfor- 
mations. Two- and three-flavor calculations of the phe- 
nomenon yield qualitatively different results. The two- 
flavor trajectory is unstable to small displacements in 
the three-flavor space. As we saw, even a tiny nonzero 
value of ArriQ is enough to give a different Ve spectrum. 
Moreover, the adiabaticity of the Amg-driven evolution 
is marginal, leading to a mixed spectrum in the chan- 
nel, rather than the usual complete swap. Clearly, a 3- 
flavor analysis is a must at this stage of the explosion. 

This letter is meant to open, rather than close, the is- 
sue. It is also important to consider what happens in the 



NH case. Briefly, the high-energy splits at « 20 MeV re- 
main, but mixed spectra appear below |49]. One should 
also investigate the sensitivity of the answer to the de- 
tails of the spectra and what happens in full multiangle 
calculations. Finally, numerous applications await ex- 
ploration: how the collective oscillations affect the signa- 
tures of the shock and turbulence, whether they change 
the r-process and the diffuse supernova neutrino back- 
ground, etc. 

We concur with the conclusions of i48j: the physics of 
supernova neutrino conversion continues to surprise us 
with its richness. 
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feedback. I am especially grateful to H. Duan for nu- 
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LDRD program of the Los Alamos National Laboratory. 
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We have also verified this with our simplified model of 
neutrino decoupling. 

Here 1/3 is the lower-energy state, following the standard 
notation for the IH. 

These two figures appear in Phys. Rev. Lett, as supple- 
mental online material, [52], because of space constraints. 
This order-of-magnitude estimate omits 0(1) factors. A 
more detailed analysis will be presented elsewhere. 
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FIG. 4: Final three-flavor conversion probabilities for neutrinos, shown in the mass basis. For example, P{i'i — >■ 1^2) indicates 
the probability, as a function of energy, that the neutrino originally in eigenstate transitions into eigenstate |i'2) by the end 
of the collective transformations (at 1000 km). The inverted mass hierarchy is assumed. The top panel specifles the labeling 
convention for the states. (This plot appears in Supplemental Online Material in the Phys. Rev. Lett, version, |52).) 
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FIG. 5: Same as Fig.|4] but for the antineutrinos. (This plot appears in Supplemental Online Material in the Phys. Rev. Lett, 
version, [52].) 



